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ABSTRACT  A  — 

A  simple  approach  based  on  Shoenberg's  theorem  is  described 
to  test  whether  a  set  of  border  points  of  a  simply  4-connected 
digital  picture  is  convex.  The  sequential  implementation  of  this 
method  is  linear  in  the  number  of  points;  the  parallel  algorithm 
needs  constant  time  only,  using  bitwise  parallel  Boolean  opera¬ 
tions  and  shifts  on  binary  matrices.  Suitable  modifications  of 
this  approach  can  be  used  for  decomposing  two-dimensional  objects 
into  convex  sets  and  for  filling  concavities. 
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1.  Introduction 


In  recent  years,  digital  image  processing  has  rapidly  grown 
into  a  major  discipline  with  applications  in  a  variety  of  areas 
[5,6].  One  important  problem  in  digital  image  processing  is  to 
test  whether  a  digital  object  is  convex  and  if  not  to  detect  the 
concavities.  There  are  a  number  of  approaches  to  this  problem 
[1,2,4,7,9,11-15].  Our  aim  here  is  to  describe  a  very  simple 
algorithm  to  test  whether  a  simply  4-connected  digital  object  S 
(i.e.,  an  object  without  holes)  has  a  border  which  is  convex  of 
not.  For  this  purpose  we  think  of  S  as  a  set  of  lattice  points 
in  the  Euclidean  plane  and  regard  it  as  a  polygon,  as  in  [8], 
with  the  border  points  of  S  as  its  vertices  and  having  edges  of 
length  1  (between  the  border  points  which  are  horizontally  or 
vertically  adjacent)  and  /2  (between  the  border  points  which  are 
diagonally  adjacent) . 

The  test  we  use  is  based  upon  a  very  powerful  theorem  due  to 
Shoenberg  [10]  which  gives  necessary  and  sufficient  conditions 
for  a  closed  polygon  to  be  convex  in  an  even-dimensional  Euclidean 
space.  In  sequential  implementation  this  approach  leads  to  a 
border  following  algorithm  which  is  linear  in  the  number  of  border 
points.  In  parallel  implementation  the  approach  leads  to  a  pat¬ 
tern  matching  algorithm,  i.e.  a  local  operation  on  binary  images, 
which  runs  within  constant  time  using  bitwise  parallel  Boolean 
operations  and  shifts  on  binary  images  as  basic  instructions. 


2 .  Shoenberg ' s  Theorem 


Let  n  =  P0,Plf...,Pk  where  Pi  =  (xii'xj.2' *  *  *  ^im1  is  a 
polygon  with  vertices  P^  (0*i*k)  in  m-dimensional  Euclidean 
space,  Em.  We  assume  that  n  spans  Em,  which  is  the  case  pro¬ 
vided  the  matrix 


X 


I  I1,xil'*i2' 

(i  =  0,1,2,. 


,k,  k=an) 


(1) 


is  of  rank  (m+1) . 

We  say  that  the  polygon  II  is  convex  on  Em  provided  it  spans 
Em  and  crosses  no  hyperplane  more  than  m  times. 

Shoenberg  [9]  proves  that  the  polygon  n  is  convex  on  Em  iff 
the  matrix  X  defined  by  (1)  is  of  rank  (m+1)  and  all  its  non¬ 
vanishing  minors  of  order  (m+1)  are  of  the  same  sign. 

This  theorem  is  intuitively  reasonable,  as  it  says  that  the 
convexity  of  n  on  E2  requires  that  no  two  triangles  P^PBPv 
(a<6<y)  have  opposite  orientations. 

The  parity  of  m,  however,  plays  a  role.  If  m  is  even,  and 
II  is  convex  in  Em,  then  also  the  closed  polygon  (with  k+1  sides) 


H1  "  P0P1 


*  PkP0 


is  convex  on  E  .  For  the  case  where  m  is  odd  it  can  be  shown 
m 

that  the  first  and  last  vector  of  a  convex  polygon  can  never 
coincide.  Fortunately,  for  our  application  in  two  dimensions 
this  theorem  is  very  useful. 


For  the  closed  polygon  H^,  it  does  not  matter  which  vertex 
is  taken  to  be  the  first  as  long  as  the  correct  cyclic  order  of 
the  vertices  is  preserved.  This  is  evident  also  from  the  fact 
that  the  cyclic  permutations  of  the  rows  of  the  matrix  X  will 
not  change  the  common  sign  of  its  minors  of  the  odd  order  of 
(m+1) .  However,  for  the  sequential  implementation  of  our  ap¬ 
proach  for  detecting  concavities  we  start  with  the  border  point 
which  is  the  first  from  left  to  right  in  the  uppermost  row  of 
the  simply  connected  digital  object  S. 


3 .  Test  for  Convexity 

The  condition  of  Shoenberg ' s  theorem  is  very  strong  (both 
necessary  and  sufficient)  to  test  for  convexity  of  polygons. 

For  two-dimensional  polygons  with  (k+1)  sides,  drawn  on  a  lat¬ 
tice  (as  in  the  case  of  a  digital  picture) ,  this  test  takes  a 
very  simple  form  of  evaluating  k+1  inner  products  of  two  three- 
component  vectors. 

Let  the  border  points  of  S  forming  a  closed  polygon  with 
(k+1)  vertices  (PqjP^, . . . ,Pk,PQ)  be  denoted  by  their  position 
coordinates  in  the  lattice,  thus: 

p0  =  *x0'y0*;  P1  =  <*i,Yi>;---;  pi  =  ^xi'yi* ; *  * ' ;Pk  =  **k'yk^ 
Then  according  to  Shoenberg’ s  theorem,  II  =  PQ...PkP0  is  a  closed 

convex  polygon  iff  the  matrix 


X  = 


1  *o  yo 


1  *i  *i 


-1  yk_J 


(2) 


is  of  rank  3  and  all  its  non-vanishing  minors  are  of  the  same 
sign. 

This  test  can  be  carried  out  as  follows: 


Let 


A.  =  Det 

l 


1  x. 


i-1  Yi-1 


1  x. 


1  xi+l  yi+l 


(3) 


_ _ : _ 


(i  is  taken  in  cyclic  order  mod(k+l),  0*ia«k)  .  Then  for  each 
i  =  0,1,..., k,  form  the  numerical  sequence  {A1# A2, . . . , Ak, AQ}  . 

If  all  of  them  are  of  the  same  sign,  II  is  convex.  Otherwise, 
there  will  be  changes  of  sign  (ignoring  zeros)  which  indicate 
the  concavities.  A  zero  value  of  A^  denotes  that  the  given 
three  points  are  collinear  (a  degenerate  triangle) . 

In  Figures  1  and  2  we  give  an  illustration  of  this  procedure. 
In  our  figures,  ■  denotes  an  object  point,  and  □  a  background 
point. 

However,  it  is  not  necessary  to  perform  the  inner  products 
of  these  three-component  vectors,  i.e.,  to  compute  the  A^'s, 
during  such  a  test  procedure.  There  exists  only  a  finite  number 
of  possible  pairs  of  such  vectors.  All  possible  combinations 
with  the  corresponding  A^-values  are  included  in  the  following 


(and  all  90°  rotations  of  these  schemes) .  In  the  first  step 
we  go  from  ^  to  the  center  point  of  such  a  scheme;  in  the  second 
step  we  go  from  the  center  point  in  a  direction  which  is  labeled 
with  value  6.  During  both  steps  the  interior  of  the  object  S 
(of  the  polygon  II)  is  on  our  right  side.  Let  the  center  point 


object  S 


polygon  II 


(AlfA2,A3,A4#A5#A0)  =  (-1,-1, 0,-1, -1,-2) 
The  object  S  (polygon  II)  is  convex. 


Figure  1.  Example  of  a  convex  object. 


object  S  polygon  II 

( A ^ , A ^ , A g , Ag , ,  Ag  )  =  0 ,  — 1, 0 , — 1, 0 , — 1, 2 , — 1) 

The  object  S  (polygon  H)  is  concave;  there  are  two  sign 
changes  in  the  cyclic  order  indicating  one  concavity  at 
point  P^. 


Figure  2.  Example  of  a  concave  object. 


Then  A. 
x 


be  the  border  point  P^. 


=6.  In  particular,  all 


possible  A^-values  are  taken  from  the  set  {-2, -1, 0 , 1, 2} . 
For  example  in  the  cases 


4 .  Sequential  Implementation 

Let  us  assume  that  in  one  binary  image  there  is  exactly 
one  simply  4-connected  object  S.  At  first,  starting  with  the 
uppermost  row  of  this  image,  we  search  down  row  by  row  from  left 
to  right  for  the  first  object  point  PQ,  which  is  necessarily  a 
border  point.  Then,  following  the  border  with  the  interior  on 
the  right  side,  we  label  all  border  points  P^,  P2'  *  *  *  ,Pk'  P0  steP 
by  step  with  A^'s  according  to  the  schemes  given  in  Section  3. 
During  this  procedure  we  can  count  the  number  of  sign  changes 
(if  there  is  any  positive  A^-value  then  the  object  is  concave) , 
or  the  number  of  border  points,  or  we  can  add  all  positive  A^- 
values  in  one  register,  say  R+,  and  all  negative  A^-values  in 
another  register,  say  R~.  At  the  end  of  this  border  following 
procedure  it  is  possible  to  use  such  results  to  compute  measures 
of  convexity — for  example 

c  ,g>  _  number  of  sign  changes  x  100 
i '  number  of  border  points 

or 

C2(S)  =  (.f+1)xl0° - 

(R++l)x  number  of  border  points 

In  Figure  3  we  give  four  examples  of  objects  and  their  measures 
C1  and  C2.  Of  course,  an  object  S  is  convex  iff  C^S)  =  0. 

Let  us  assume  the  object  S  is  encoded  in  Freeman  code  repre¬ 
sentation  (see,  e.g.,  [4])  where  the  directions  are  encoded 
according  to  the  following  scheme: 


0 


c2(Sl) 


=  112.5000 


s 


2 


Figure  3.  Examples  for  C-^  and  C2  values 


For  example,  the  Freeman  code  representation  for  object  S4  in 
Figure  3  is  0002334765.  It  is  possible  to  calculate  the  A^- 
sequence  for  a  given  object  S  (and  related  values,  like  R+  or 
R  )  while  reading  the  Freeman  code  sequence  of  the  object  from 
left  to  right,  using  the  labeling  schemes  given  in  Section  3. 


For  example,  for  the  object  S4  we  get 

00  -*■  Ax=0  23  f  A4=-l  47  -*■  A?=-l  50  A  =-l 

00  A2=0  33  -*■  A5=0  76  Ag  =  l 

65  *  A_=l 


02  -  A,=-l 


34  -*■  Ar=-1 


5 .  Parallel  Implementation 


Let  us  assume  that  in  one  binary  image  of  size  NxN  there  is 
exactly  one  simply  4-connected  object  S.  We  choose  as  a  irodel 
for  parallel  computation  the  PBS  parallel  binary  image  proces¬ 
sing  system) ;  the  exact  definition  of  this  model  can  be  found  in 
[3] .  The  PBS  can  perform  bitwise  parallel  Boolean  operations  and 
shift  operations  on  binary  images  as  basic  instructions.  For 
example,  the  program  for  computing  the  edge  image  for  the  given 
object  S  has  the  following  form: 

input  binary  image  X  of  size  NxN  with  object  S; 
begin  binary  images  Y,Z  of  size  NxN; 

Y=X-*-l;  Z=XAY;  Y=X«-1;  Z=ZAY; 

Y=Xtl;  Z=Z  AY ;  Y=Y  +  1;  Z=ZAY; 

Z=Z;  Y=XAZ;  print  Y 

end. 

In  [3]  it  was  explained  that  the  PBS  is  very  time  efficient  for 
computing  local  operations  on  binary  images;  edge  detection  is 
one  example  of  such  a  local  operation.  Fortunately,  the  label¬ 
ing  of  the  border  points  of  S  with  A^-values  is  such  a  local 
operation. 

Essentially,  the  computation  of  the  A^-values  can  be  described 
as  a  special  pattern  matching  procedure.  The  necessary  size  of 
the  patterns  is  3x3  only.  In  Figure  4  we  have  listed  all 
patterns  which  must  be  matched  during  this  procedure.  In  this 
figure,  ■  represents  an  object  point,  □  represents  a  background 


>w< 


an  object  point  or  a  back- 


point,  and  Ei  can  be  either,  i.e., 

ground  point.  If  a  pattern  matches  a  3x3  window  in  the  binary 

image,  the  center  point  of  this  window  is  labeled  with  the 

1 

corresponding  A^-value.  Using  binary  images  in  this  procedure 
only,  this  labeling  can  be  done  in  several  bit  planes,  one,  for 
example,  corresponding  to  label  1,  and  one  to  label  -1.  This 
labeling  can  be  done  on  PBS  within  constant  time.  After  this 

procedure,  we  can  count  the  l's  in  the  resulting  bit  planes 

2  +  — 
within  time  O((log  N)  )  (see  [3])  and  get  the  values  R  and  R 

defined  in  Section  4,  for  example. 

Let  us  assume  the  object  S  is  encoded  in  Freeman  code  repre¬ 
sentation  already,  igi^^. .  .i^.  Assume  that  k+1  processors  are 
available  for  parallel  computation.  In  a  first  step,  each  pro¬ 
cessor  with  number  j  (0* j^k)  gets  the  values  i^,ij+^,  as  input 
and  computes  A^.  In  a  second  step  each  processor  compares  its 
value  with  the  values  of  its  two  neighbor  processors.  Thus,  after 
two  steps  we  discover  whether  the  object  is  convex  or  not.  The 
values  R+  and  R-,  or  the  number  of  sign  changes,  can  be  computed 
using  k+1  processors  within  llog2kJ  steps. 


1In  using  these  masks,  we  assume  that  the  object  thickness  is 
at  least  two;  otherwise,  pattern  matching  with  these  masks  may 
result  in  non-unique  labels. 


-v - Nfe. 


*  * 


6 .  Concluding  Remarks 

a.  The  approach  described  for  detection  of  concavities  can 
be  applied  to  decompose  polygons  into  convex  sets. 

b.  It  can  also  be  used  to  fill  concavities  by  addition  or 
deletion  of  edges  at  the  appropriate  vertices. 

c.  An  iterative  application  (relaxation  algorithm)  of  com¬ 
parisons  between  A ^-values  in  the  neighborhoods  of  border  points 
can  be  used  for  smoothing  of  "small  concavities." 

d.  After  this  report  was  prepared,  the  authors  came  across 
a  recent  paper  by  Bribiesca  and  Guzman  [16]  in  which  a  method 
based  on  Freeman  Code  is  given  for  the  description  as  well  as 

measurement  of  differences  in  shapes  (shape  taxonomy)  among  poly¬ 
gonal  boundaries — a  shape  being  defined  as  that  property  of  an 
object  which  is  invariant  under  translation,  rotation  and  similar¬ 
ity  transformation  (choice  of  grid  or  lattice  size  and  its  orien¬ 
tation)  .  In  order  to  arrive  at  a  canonical  description,  the 
orientation  of  the  grid  is  normalized  by  taking  one  of  its  axes 
parallel  to  the  length  of  a  basic  rectangle  which  just  contains 
the  boundary;  however,  the  choice  of  grid  size  still  plays  a  role 
in  uniqueness.  For  this  purpose  the  following  procedure  is  used; 
each  lattice  point  along  the  boundary  is  encoded  in  ternary  base, 
viz .  0  for  Convex,  1  for  Straight,  and  2  for  Concave.  This 

resulting  ternary  number  is  then  circular  shifted  until  it  assumes 
the  least  value;  this  serves  as  the  shape  number;  the  number  of 


digits  chosen  naturally  depends  on  the  grid  size  and  hence  for 
uniqueness  the  number  of  digits  (called  order)  is  prescribed. 

It  is  clear  that  the  algorithm  described  in  our  paper  can 
be  suitably  modified  if  we  use  4-connectedness  for  parallel 
automatic  generation  of  such  shape  numbers. 

e.  Limitations  of  this  approach:  Finally  a  note  of  caution: 
in  the  method  described  here,  we  assume  that  the  given  object 
is  digitized  and  the  resulting  polygonal  boundary  is  examined 
for  convexity.  It  may  happen  that  new  concavities  are  introduced 
in  digitization  (due  to  orientation  and  size  of  mesh)  which  may 
not  exist  in  the  original  object.  For  example,  the  rotation  of 
a  square  when  digitized  will  exhibit  concavities.  This  limita¬ 
tion,  however,  exists  in  other  available  methods  also  (16] . 


References 


[1]  Arcelli,  C.,  and  Cordelia,  L.,  Concavity  point  detection  by 
iterative  arrays.  Computer  Graphics  Image  Processing,  3, 

pp.  34-47,  1974. 

[2]  Hodes,  L.,  Discrete  approximation  of  continuous  convex  blobs, 
SIAM  J.  Appl.  Math.,  19,  pp.  477-485,  1970. 

[3]  Klette,  R. ,  Parallel  operations  on  binary  images.  Computer 
Science  Center,  University  of  Maryland,  College  Park,  TR-822, 
November  1979. 

[4]  Pavlidis,  T.,  Structural  Pattern  Recognition,  Springer-Verlag, 
New  York,  1976. 

[5]  Rosenfeld,  A.,  Digital  Topology,  Amer.  Math.  Monthly  86 , 
pp.  621-630,  1979. 

[6]  Rosenfeld,  A.,  and  Kak,  A.  C.,  Digital  Picture  Processing, 
Academic  Press,  New  York,  1976. 

[7]  Rutovitz,  D.,  An  algorithm  for  in-line  generation  of  a  convex 
cover.  Computer  Graphics  Image  Processing,  £,  pp.  74-78,  1975. 

[8]  Sankar,  P.  V.,  and  Krishnamurthy,  E.  V.  ,  On  the  compactness 
of  subsets  of  digital  pictures.  Computer  Graphics  Image 
Processing,  8^  136-143,  1978. 

[9]  Schachter,  B. ,  Decomposition  of  polygons  into  convex  sets, 

IEEE  Trans.  Comp.,  C-27 ,  pp.  1078-1082,  1978. 

[10]  Shoenberg,  I.  J.,  Isoperimetric  inequality  for  closed  curves 
convex  in  even  dimensional  Euclidean  space,  Acta  Mathematica, 
91,  pp.  143-164,  1954. 

[11]  Sklansky,  J.,  Chazin,  R.  L.,  and  Hauser,  B.  J. ,  Minimum 
perimeter  polygons  of  digitized  silhouettes,  IEEE  Trans. 

Comp.,  C-21 ,  pp.  260-268,  1972. 

[12]  Sklansky,  J.,  Measuring  concavity  on  a  rectangular  mosaic, 

IEEE  Trans.  Comp.,  C-21,  pp.  1355-1364,  1972. 

[13]  Sklansky,  J. ,  Recognition  of  convex  blobs.  Pattern  Recognition, 
2_,  pp.  3-10,  1970. 


[14]  Sklansky,  J.,  On  filling  cellular  concavities.  Computer 
Graphics  Image  Processing,  £,  pp.  236-247,  1975. 

[15]  Sklansky,  J.,  Cordelia,  L.  P.,  and  Levialdi,  S.,  Parallel 
detection  of  concavities  in  cellular  blobs,  IEEE  Trans. 

Comp.,  C-25 ,  pp.  187-196,  1976. 

[16]  Bribiesca,  E.  and  Guzman,  A.,  How  to  describe  pure  form 

and  how  to  measure  differences  in  shapes  using  shape  numbers, 
Proc.  IEEE  Conf.  on  Pattern  Recognition  and  Image  Processing, 
Chicago,  Illinois,  pp.  427-436,  August  1979. 


SECURITY  CLASSIFICA 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS  . 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO.  3.  RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (and  Subtitle) 


BBLSiW 


ALGORITHMS  FOR  TESTING  CONVEXITY  OF  DIGITAL 
POLYGONS 


5.  TYPE  OF  REPORT  &  PERIOD  COVEREO 


Interim 


6.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR!®.) 

Reinhard  Klette 
E.  V.  Kr ishnamurthy 


9.  PERFORMING  ORGANIZATION  NAME  ANIVADDRESS 

Computer  Vision  Laboratory,  Computer- 
Science  Center,  University  of  Maryland 
College  Park,  MD  20742 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


Air  Force  Office  of  Scientific  Research/NM 
Bolling  AFB,  Washington,  D.  C.  20332 


14.  MONITORING  AGENCY  NAME  a  ADDRESS!//  dillerent  from  Controlling  Office)  I  15.  SECURITY  CL  ASS.  (ot  this  report) 


8.  CONTRACT  OR  GRANT  NUMBER!®.) 


AFOSR  77-3271 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  a  WORK  UNIT  NUMBERS 


61102F  2304/A2 


12.  REPORT  DATE 

December,  1979 


13.  NUMBER  OF  PAGES 
20 


UNCLASSIFIED 


ISa.  DECL  ASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  this  Report) 

Approved  for  public  release;  distribution  unlimited. 


>7.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 


19.  KEY  WORDS  ( Continue  on  reverse  side  If  necessary  and  Identify  by  block  number) 

Pattern  recogintion 
Image  processing 
Shape 
Convexity 
,  Polygons 
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A  simple  approach  based  on  Shoenberg's  theorem  is  described  to  test  whether  a 
set  of  border  points  of  a  simply  4-connected  digital  picture  is  convex.  The 
sequential  implementation  of  this  method  is  linear  in  the  number  of  points;  the 
parallel  algorithm  needs  constant  time  only,  using  bitwise  parallel  Boolean 
operations  and  shifts  on  binary  matrices.  Suitable  modifications  of  this 
approach  can  be  used  for  decomposing  two-dimensional  objects  into  convex  sets 
and  for  filling  concavities. 
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